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Abstract: We continue our investigations on the relation between hydrodynamic and 
higher quasinormal modes in the AdS black hole background started in arXiv:0710.4458 
[hep-th]. As is well known, the quasinormal modes can be interpreted as the poles of 
the retarded Green functions of the dual N = 4 gauge theory at finite temperature. The 
response to a generic perturbation is determined by the residues of the poles. We compute 
these residues numerically for energy-momentum and R-charge correlators. We find that 
the diffusion modes behave in a similar way: at small wavelengths the residues go over 
into a form of a damped oscillation and therefore these modes decouple at short distances. 
The sound mode behaves differently: its residue does not decay and at short wavelengths 
this mode behaves as the higher quasinormal modes. Applications of our findings include 
the definition of hydrodynamic length and time scales. We also show that the quasinormal 
modes, including the hydrodynamic diffusion modes, obey causality. 
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1. Introduction 

In recent years a new paradigm concerning the high temperature behaviour of QCD has 
been established: the strongly coupled Quark-Gluon Plasma (sQGP). Experimental re- 
sults from heavy ion collisions at RHIC indicate that QCD at temperatures around 2T C 
is strongly interacting, in spite of being in a deconfined phase, and thus rendering per- 
turbative computations not suitable for describing it. While static properties of strongly 
coupled gauge theories at finite temperature can be readily analyzed on the lattice, the 
study of out-of-equilibrium phenomena faces considerable difficulties. In the last years the 
AdS/CFT correspondence [1-3] has emerged as a useful tool to understand analytically 
the dynamics of non-Abelian gauge theories in the strongly coupled plasma phase. 

According to AdS/CFT, the asymptotically AdS black hole is dual to the plasma phase 
of the strongly coupled J\f = 4 gauge theory [4-6] . The line element is 

d S 2 = ^(-/(r)d^ + dx^)+^, (1.1) 
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with f(r) = (1 — r 4 /r 4 ). The Hawking temperature Tjj = tq/ttL 2 is the temperature T 
of the dual field theory. Real-time correlation functions in the thermal gauge theory can 
be computed using classical solutions of fields living in the gravity dual [7]. If a field is 
excited in the presence of a black hole, the energy of the fluctuation will be lost inside 
the horizon and eventually the final state will be a larger black hole with no fluctuations. 
This process is described by the quasinormal spectrum, that was first computed for black 
holes in asymptotically Anti de Sitter (AdS) spacetimes in [8]. Quasinormal modes (at 
fixed momentum) exist only for a discrete set of complex frequencies that can be identified 
with the poles of the retarded Green functions or resonances in the dual field theory [9-11]. 
Therefore, quasinormal modes describe dissipation processes in the plasma. 

The retarded Green functions depend on frequency uj and momentum q. The quasi- 
normal modes can be understood as solutions of (Gr) -1 = in the complexified u> plane. 
However, one can also ask for solutions in the complexified momentum plane. Again one 
finds a discrete spectrum of complex momentum modes (keeping the frequency real) whose 
imaginary part can be interpreted as the inverse absorption length [12]. We will also 
consider these complex momentum modes and their residues in this paper. 

Hydrodynamic modes like diffusion, shear or sound modes are also described by quasi- 
normal modes in the gravity dual. They are special cases of quasinormal modes whose 
frequencies vanish in the zero momentum limit [13, 14]. One of the most interesting results 
of this AdS black hole hydrodynamics has been the derivation of a universal bound for the 
shear viscosity to entropy ratio 2 > ^JL_ [15 ; 16]. It has also been argued that this value 
is relevant to the description of heavy ion phenomenology at RHIC [17]. The shear vis- 
cosity is a hydrodynamic transport coefficient governing the momentum diffusion through 
the medium. One way to derive it is to compute the lowest quasinormal frequency in the 
retarded two-point correlator in the vector channel of the stress tensor. 

The hydrodynamic approximation can be interpreted as a low energy effective theory 
where all the degrees of freedom except the hydrodynamic modes have been integrated out. 
It is then possible to make a systematic expansion in higher derivative terms, whose coef- 
ficients depend on the microscopic physics. The first order terms correspond to the usual 
hydrodynamic equations, while second order terms take into account the finite delay in 
the response of the system. Second order hydrodynamics is important since the relativistic 
first order formalism leads to the well-known problem of acausal signal propagation. As we 
will review in appendix [C], the speed of signal propagation is determined by the so called 
front velocity v = lim^^oo with q R the real part of the complex momentum mode. 
Diffusion equations as they arise typically in first order hydrodynamics do not have a finite 
front velocity and therefore allow acausal behaviour. In second order hydrodynamics the 
non-zero delay of the response is encoded in a relaxation time constant and can restore 
causality [18-21]. The latter have been recently computed for conformal theories using 
the AdS/CFT correspondence [22-25]. However, even second order hydrodynamics cannot 
simply be extended to arbitrary short wavelengths or high frequencies and thus a front 
velocity smaller than the speed of light in second order hydrodynamics is not guaranteed, 
since inevitably the hydrodynamic description breaks down at small wavelengths. Thus 
the question arises how causal signal propagation is guaranteed in the N = 4 plasma at 
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strong coupling. We will show numerically that the hydro dynamic mode approaches the 
front velocity vp = 1 and argue that this is indeed the exact value for all quasinormal 
modes including the hydrodynamic diffusion modes. 

As an effective theory, the hydrodynamic approximation will be valid up to energies 
of the order of the microscopic scale, in this case given by the temperature T. At weak 
coupling hydrodynamic modes are collective excitations made up by the collective motion 
of the hard particles. Therefore they actually should decouple at short wavelengths. In 
standard perturbative finite temperature field theory and in the HTL approximation col- 
lective excitations do indeed appear already in the two-point correlators of the fundamental 
fermion and gauge fields. At weak coupling one distinguishes a soft scale of order gT from 
a hard scale T. Above T the hard partons are the quarks and gluons of the elementary 
fields, at the soft scale there are quasiparticles, dressed quarks and gluons and collective 
excitations, i.e. the plasmon and plasmino modes [26-28]. They show up as poles in the 
retarded two-point functions but can be distinguished by the behaviour of the residues at 
momenta at or above the hard scale T: the residues of the collective excitations decay 
exponentially signalling their decoupling at short wavelengths. What can we expect then 
for the quasinormal modes? In fact the high frequency and short wavelength limit can be 
interchanged with the T — > limit [10]. Going to zero temperature an infinite number of 
quasinormal modes come together and open up a branch cut, characteristic for the T = 
retarded correlators. Therefore, the residues of almost all quasinormal modes should not 
decay in this limit. However, a finite number might still decouple, and indeed this is what 
we find for the residues of quasinormal modes representing diffusive behaviour at long 
wavelengths. 

In this work we extend previous analysis of the residues of the quasinormal modes of 
the R-current [29] to energy-momentum tensor correlators. In section [2] we recall linear 
response theory and how it is related to hydrodynamic behaviour. We study the regime of 
validity of the hydrodynamic approximation and define a lower bound in length and time 
scales. In section || we discuss how the retarded Green function can be expressed as infinite 
sums over the poles at the quasinormal frequencies or complex momenta. We point out 
that these sums do not converge and analytic terms are needed to regularize the sums. 
These analytic terms should not be confused with contact terms arising from the action. 

In section |] we recall and extend the results of ref. [29] on R-charge correlators. In 
particular we compute the residues of the complex momentum modes whose imaginary 
parts give the inverse absorption lengths [12] and study how well can the real and imaginary 
parts of the retarded Green functions be approximated by the contributions of the lowest 
quasinormal modes. 

In section [| we study quasinormal frequencies, complex momentum modes and their 
residues for stress tensor correlators. We find that the residues of the shear mode quasi- 
normal frequency go over into a form similar to a damped oscillation at short wavelengths. 
Numerically we find that it decouples for short wavelengths. The behaviour of the quasi- 
normal frequencies in the shear channel are very similar to the diffusion channel for the 
R-charge correlators. The hydrodynamic mode crosses the imaginary parts of the other 
(low) quasinormal frequencies roughly at the locations of the zeroes in the shear mode 



- 3 - 



residue, signalling a breakdown of the effective field theory based on the shear mode alone. 
We also find that the residues in the sound channel show a different behaviour: the lowest 
mode, the sound mode, does not decouple at short wavelengths. It becomes similar to a 
higher quasinormal mode. Again we study how well the spectral functions can be approx- 
imated by summing only over the lowest quasinormal modes. In addition we also compute 
the residues in the shear and sound channels of the complex momentum eigenvalues for 
fixed real frequencies. This is important in order to compute (numerically) the front veloc- 
ity. We find that the front velocity in all modes, even the hydrodynamic one, approaches 
unity at high frequency and give an argument that this is indeed enforced by the Lorentz 
symmetry of the underlying T = theory. 

In appendix [A] we collect technical details on how we compute the retarded Green 
functions, their poles and residues. In appendix [B| we show how the location of the zeroes 
in the residues of R-charge diffusion and shear mode are determined by algebraic equa- 
tions arising from the recursion relations of the corresponding Heun equations. Finally in 
appendix ^ we recall the definition of the front velocity. 

2. Hydrodynamic scales and linear response theory 

We consider a medium in thermal equilibrium. To study the effect of small external per- 
turbations it is enough to use linear response theory when the energy of the perturbation 
is negligible compared to the total energy of the system. In this linear approximation, the 
response of a field $ to a perturbation represented by an external source j(t,x) is 

(*(t,x)> = - JdTd 3 £G a (t-T,x-t)j(T,t), (2.1) 
where Gr is the retarded two-point correlation function, 

G R (t - r, x - = -i 9(t - r) <[$(t, x), <D(r, £)]> • (2-2) 
Using the Fourier transforms of the Green function and the source this is rewritten as 

<*(*,*)) = - j e-- f +^G R (u,,q)j>,q) . (2.3) 

We can now make the analytical continuation to the complex uj plane and use Cauchy's 
theorem to evaluate this integral by closing the contour on the lower half of the complex u 
plane for t > 0, thus picking up the contributions from the poles of the retarded propagator, 
i.e. the frequencies of the quasinormal modes in the AdS/CFT correspondence. At this 
point we assume that the retarded Green function is analytic in the upper half of the 
complex frequency plane and that its only singularities are single poles in the lower half 
plane. This is indeed the case for the holographic retarded two-point functions. We also 
assume that the source does not introduce new non-analyticities in the lower half plane. 
The retarded propagator is therefore of the form 

G R ( W , q )^^^, (2.4) 

poles 
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where io n := Q n — iT n and R n are the residues of the retarded propagator evaluated at the 
poles. We should emphasize that this is a formal expression since the infinite sum does not 
necessarily converge; we will come back to this issue in section ||[ Doing a partial Fourier 
transform in space, the response of the system to a perturbation is now expressed by a sum 
over the poles 

<<&(t,q)> = i^(t)^i? n (a; n ( (? ),q)JK(q),q)e-^^- r "^ . (2.5) 

n 

The location of the quasinormal modes determine their frequency and damping and their 
residues determine how much each mode contributes to the response. The quasinormal 
mode spectrum of the planar AdS black hole is well known by now [10,30-32]. For each 
channel corresponding to a gauge invariant operator in the dual gauge theory there is an 
infinite tower of quasinormal modes. Since for holographic duals of gauge theories there 
exist infinitely many quasinormal modes the precise form of the response also depends on 
the form of the source. There is however an exception to this: the extreme long-wavelength 
modes of conserved charges such as energy and momentum. Perturbations that induce a 
change in the overall charge have to excite the special hydrodynamic modes. These are 
modes whose quasinormal frequencies obey lim q _>o wh (q) = 0. This behaviour guarantees 
that the integrated response f d 3 x {$>(t, x)) is time independent, as it has to be for a 
conserved charge. Notice however that there is no fundamental reason that would forbid 
neutral perturbations to dissipate away primarily in the higher quasinormal modes, i.e. 
the infinite tower of modes with non vanishing zero momentum limit. It is also interesting 
to think about what the presence of the infinitely many quasinormal modes mean for the 
formulation of an initial value problem. Let us assume for the moment that there would 
be only one quasinormal mode, e.g. a diffusion mode in a conserved charge. The time 
development for t > would then be completely determined by specifying the expectation 
value of the field at t = 0, i.e. (<I>(0,x)) = <£o( x )- Indeed, cutting the sum in ( |2.5| ) after 
the first term, the hydrodynamic term in this case, leaves us with a one-to-one relationship 
between the initial value and the source $o(ci) = -R(q)j(q)- However, in our case there 
are infinitely many quasinormal modes and therefore one has to specify not only the field 
at t = but also an infinite number of its time derivatives in order to be able to compute 
the source from the initial condition. Strictly speaking the system becomes non-Markovian 
due to the presence of the infinitely many quasinormal modes: it remembers (at least for 
short times) the history of how it reached a certain state. Putting it another way, we can 
define a system to be in local thermal equilibrium if it has no memory, i.e. the response is 
completely dominated by the hydrodynamic mode. 

We will concentrate on simple forms of source perturbations. We imagine that the 
source acts only over a time interval At. Such a perturbation will excite a significant 
number of quasinormal modes for small At and in the limit where At — > it will excite 
all quasinormal modes with equal weight. We can then assume a source of the form 
i(t,x) = £(i)cos(qx). 

At small values of the momentum, the hydrodynamic mode dominates the long time 
behaviour since it has the smallest imaginary part. In fact, we can define the hydrody- 
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namic time scale, from which on the hydrodynamic mode dominates and the hydrodynamic 
approximation is good, by demanding that the response in the hydrodynamic mode equals 
the response in the first quasinormal mode at the time rjj. We can estimate this time scale 
using (2.5) 

log |i?H Jh| - log|i?i h\ 

th = r— ^ , (2.6) 

where we have written j n = J(w n (q),q). The hydrodynamic description will be trustable 
for times larger than th , whereas for shorter times the contribution of higher modes must 
be taken into account. 

Now we switch the roles of time and one spatial coordinate and choose a periodic 
perturbation localized in space, over an interval Ax in the x direction. In the limit Ax — ► 
we assume a source of the form j(t,x) = 8(x)exp[— i(ut — kj_x_j_)]. Prom ( |2.1[ ), the response 
of the system takes the form 

($(i,x)) = -i-e-^'- k ^ J dge^GRfok^g) . (2.7) 

Such a perturbation has the form of a plane wave in the perpendicular directions xj_. We 
have assumed that the perturbation started far in the past such that the system has reached 
a stationary state. We will also assume that it is no further modulated in the xj_-directions, 
i.e. we set kj_ = 0. By Cauchy's theorem we can evaluate this integral by closing the contour 
on the upper half of the complex q plane for x > (on the lower-half complex q plane for 
x < 0). The contributions from the poles of the retarded propagator in the complex q 
plane correspond to complex wavenumber modes in the AdS/CFT correspondence, whose 
imaginary part determine the absorption lengths of perturbations in the plasma [12,33]. 
We consider again that the only singularities of the retarded propagator are single poles in 
the complex momentum plane, so 

,9- <?nM 
poles 

where q n = q^ + iq^. The response of the system to a perturbation localized in space as a 
sum over the poles reads 

($(t, x)> = -i sign(x) e~ lujt K( u > e iq * (u;)x - qljuj)x . (2.9) 

n 

By parity symmetry x — > —x, if there is a pole at q n = q^ + Wm then also q = —q n has to 
be a pole, so the poles of the retarded Green function lie in the first and third quadrants 
of the complex q plane for x > and x < 0, respectively. On the supergravity side this has 
been indeed proved in [29]. 

Analogously to what we did before, we can define a hydrodynamic length scale. From 
Q we find 

_ WI-MBl! (2 . 10) 

Im qn — Im q\ 
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For distances from the origin of the perturbation larger than in, the hydrodynamic mode 
will dominate the response of the system. We emphasize that the definition of the hy- 
drodynamic length scale ( 2.10| ) applies only if the sources are localized in space. For a 



more general source, the Fourier transform evaluated at the complex momentum mode 
j(uj , q n (uj)) will enter in the definition of £u in an analogous way to ( |2.6| ). 

As we have seen the response and the hydrodynamic scales depend crucially on the 
knowledge of the residues of the retarded correlators. In the following we will determine 
them numerically and study such hydrodynamic scales. 

3. Including higher thermal resonances 

The hydrodynamic mode provides a good description of the system at very long wavelengths 
u, q <C T. As we have argued, for shorter times or distances the holographic computation 
shows that higher modes start to be relevant. It is an interesting question to see if a 
description in terms of the hydrodynamic mode plus a few thermal resonances gives a 
reasonable good model of the plasma up to frequencies of the order of the temperature 
lo ~ T. What we mean by thermal resonances are the quasinormal modes found in the 
holographic dual, that describe the dissipation of gauge-invariant configurations. 

Due to conformality, the properties of the N = 4 plasma rescale trivially with the 
temperature T. A possible way to introduce a non-trivial temperature dependence would 
be to compactify the theory in a three sphere of radius R. The compactification breaks 
conformal invariance, so there is a non-trivial dependence on RT. Then, the physics of the 
plasma in flat space R — > oo can also be recovered in the infinite temperature limit T — > oo 
[6]. Computations of the quasinormal mode spectrum at lower temperatures RT > 1 
[8, 34] show that the lowest modes are long lived enough to have a good quasiparticle 
interpretation. Other computations in the context of flavour branes, where high and low 
temperature is given in terms of the quark mass m q /T, also show similar results for the 
quasinormal spectrum of mesons [33,35-37]. This suggests that for more realistic plasmas, 
e.g. non-conformal and maybe closer to the intermediate temperature regime between the 
deconfinement transition and the free gas limit than the simpler models just mentioned, 
a description in terms of bound states could be appropriate, although usually coloured 
bound states are considered [38-40]. In spirit, this approach is similar to the description 
of QCD correlators at low energies using the lowest states of the meson spectrum. 

From the linearized AdS/CFT computation, we know that the retarded correlators of 
channels with hydrodynamic modes have the general form 

where Oh, w n are the hydrodynamic and quasinormal poles, Rn, R n their residues and 
«h is the number of hydrodynamic modes. For the shear and diffusive channels «h = 1, 

(2) (1) * (2) (1) * i 

while for the sound channel nn = 2 and f% = — , i?y = — R^ . The term A 



1 There exist of course also channels without any hydrodynamic behaviour nn = 0, such as scalar field 
perturbations. We will not consider these channels in this paper. 
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denotes possible terms analytic in frequency and momentum. In principle they look like 
the contact terms that can arise from Ward identities or some choice of renormalization 
scheme, in which case they could be absorbed by a finite number of counterterms. However, 
their origin is different and we will see that they are necessary in order to give a consistent 
Green function, the main reason being the infinite sum over quasinormal modes. 

Let us concentrate on a particular example, the density-density correlator of the re- 
current. In [36] the explicit form of the current-current correlator at zero momentum was 
found (to := w/(2irT)), 



N 2 T 2 { 
Gij(to) = s ij—g — \ ixv + to 2 



. (1-i) \ , f-(l + i) 



^ to - to n 



The poles and its residues can be extracted from the following representation of the 
digamma function 



Y>(*) = -7E - E ( \— - -) , (3.3) 

71=1 



so the residues will scale as R n ~ ro 2 r n , with r n ~ 1 and the frequencies ro n ~ n. Evaluated 
at the pole this gives R n ~ ro 2 r n ~ n 2 , as was confirmed numerically in [29]. Using current 
conservation, the density-density correlator is related to the longitudinal current-current 
correlator through Gu = jjzG 1 ^, so the residues of the density-density correlator go as 
Rn ~ q 2 r n ~ q 2 - 

The imaginary part of the retarded correlator receives a contribution from the quasi- 
normal modes 

^(^-Re^) 2 + (Im^) 2 ^n 2 

The sum over the quasinormal modes is convergent in this case. On the other hand the 
real part goes as 

Re G R ~ y K^RgKjn) 2yl (3 _ 5) 



n 



In this case the sum is not convergent, although the Green function is finite. The divergence 
of the sum over quasinormal modes can be cured by subtracting a divergent analytic part, 
or in a better defined way, subtracting order by order a small analytic part. The analytic 
term would look 

A = yA n ~ V£i- (3.6) 



n 

n n 



In order to give a consistent approximation to the retarded correlator, the analytic term 
must be properly taken into account. For the first m modes 



Gr(uj, q) ~ — . 2 , - + q > p-r H ■ — —-^ - C m q , u,q < 1 , (3.7) 



n=l 
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where we have defined the hydrodynamic pole and residue in terms of a g-dependent 
'conductivity' and 'diffusion coefficient' Ru = —iq 2 <T}i(q) and fin = —iq 2 Dn(q)- In the 
zero momentum limit, they take the value of the known constant transport coefficients 
D = l/(27rT), a = N 2 T/16ir. The coefficient of the analytic term C m depends on the 
number of modes and diverges when m —* oo. The exact Green function including all the 
modes would differ by a small amount from this approximation at small enough frequencies. 

In the hydrodynamic approximation only the mode associated to diffusion would be 
considered, adding higher derivative corrections to the effective description as momentum 
increases. It is then interesting to find out its exact behaviour beyond the hydrodynamic 
limit from the holographic dual. In ref. [29], the quantities crn(q) and Dn(q) were computed 
numerically. It was found that the residue of the diffusion mode is oscillatory, so for a finite 
value of q < 2irT, the ^-dependent conductivity we have defined becomes zero and then 
negative. Notice that the analytic corrections that can come from higher modes would not 
solve this problem, so clearly the diffusion mode does not have a good hydrodynamical 
interpretation anymore. Furthermore, the residue vanishes exponentially with increasing 
momentum, so this mode decouples and the system leaves the hydrodynamic regime. In 
the section |5| we will find an analogous behaviour in the shear diffusion mode of the stress 
tensor two point function and explain its properties there in more detail. 



4. More on R-charge diffusion 

We now consider the conserved current associated to the global R-symmetry in the 
j\f = 4 theory 

d% = . (4.1) 

The diffusion of the R-charge Q = J d 3 x Jo is described by a hydrodynamic mode, that at 
low frequencies and momenta dominates the response of the system to small perturbations 

(6J„(t,x)) = - J dt'dVG^-^x-xOr^x') . (4.2) 

The global current is dual to a five-dimensional gauge field Am- 2 We work in the gauge 
A r = 0. The gauge invariant quantities are the longitudinal and transverse electric fields 
El = qAo + uj q • A/q, = wA^, where q • At = 0. We describe the method to find the 
quasinormal modes in the appendix |A[ 

The retarded two-point functions can be expressed in terms of two scalar functions TIt 
and TLl corresponding to transverse and longitudinal polarizations, assuming q = (0,0, q) 

Gtt = IIt , G t t 
G tz (uj,q) = -G tt (u>,q) , G zz (uj,q) 



2 That is actually a component of the metric with the group index associated to an internal space. 



uj 2 — q 2 



(4.3) 



-oG tt (u,q) 



(4.4) 
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The quasinormal modes and residues R n ' of the scalar functions TLt,l were computed 
in [29]. We remark that the residues of the density-density correlator are actually 

In [29], we found that the residue of the diffusion mode shows a damped oscillatory be- 
haviour with the momentum. The zeroes appear when the diffusion pole reaches the value 
of the imaginary part of the quasinormal frequencies at zero momentum. The values of 
the momentum when this happens can be found analytically (appendix [B]). At a slightly 
smaller value of the momentum the imaginary part of the n th higher quasinormal mode 
becomes smaller than the value for the hydrodynamic mode, i.e. the diffusion mode crosses 
the n th quasinormal mode. The time development is then dominated by the contribution 
of the lowest quasinormal modes instead of the hydrodynamic mode. 

A way to estimate the regime of validity of the hydrodynamic approximation is to intro- 
duce an external perturbation, localized in time or space and compute the hydrodynamic 
time ( |2.6| ) or length ( 2.10| ) scales. In the first case we consider the evolution of the charge 
density (given by Gu), while in the second case we do it for the longitudinal component of 
the current (given by G zz ). The time- localized source is a plane wave of fixed momentum 
q that begins at t = and last some time At. The perturbation creates inhomogeneities 
that start to relax to equilibrium after it is switched off. The late time relaxation will be 
dominated by the diffusion mode. If the perturbation lasts for a long time, it is possible for 
the system to reach a steady state already dominated by the diffusion mode. This seems 
to be the case as long as At is a few times the inverse temperature. If the perturbation 
is short lived, then we find that the minimal time to enter the diffusion dominated regime 
is 27tTth ~ 0.35. We also find that for short wavelengths, when q approaches the value 
of the zero of the hydro residue, rjj grows unbounded. So for wavelengths of the order of 
the inverse temperature or higher, the relaxation to equilibrium is never dominated by the 
diffusion mode, but by higher modes. 

The space-localized source is a pulse constant on x and y directions and localized in 
the z direction with size Az. The pulse is oscillating with some fixed frequency u. The 
perturbation will be screened as we move in the z direction far from the plane where 
the pulse is localized. At large enough distances, the diffusion complex momentum mode 
would dominate the decay of the perturbation. We find that if the size of the pulse is 
larger than the scale given by the temperature Az > 1/T and the frequency is low enough 
lo < 2ttT, the screening is dominated by the hydro mode at any distance. However, 
at higher frequencies the diffusion mode starts dominating at larger distances from the 
source, so £u > 0. When the size of the source is of the order of the inverse temperature, 
finite size effects like the shape of the perturbation start to be important to determine the 
hydrodynamic length scale. In our definition of hydrodynamic time and length scales there 
is an explicit dependence on the source, 

sin(uAt) i oca ji zec [ [ n time 

(4.6) 

sin(gAx) l oca lj zec [ j n S p ace 
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The source is evaluated at the value of the of the mode u; n (q) or q n (w). For large values 
of the frequency or the momentum, the ~ 1/to, 1/q scaling determines the response, while 
for very low values it is constant. For intermediate values however, the sine function 
introduces an oscillatory behaviour. In general, we will find different answers for different 
kind of perturbations, depending on their shape. 

In summary, no matter the value of q for the time-dependent source there is a minimal 
hydrodynamic time scale th if the duration of the perturbation is short enough. However, 
for the space-localized source one can always lower the frequency to a value where the 
length becomes zero, no matter how small the source is. 

The validity of the hydrodynamic regime can also be studied through the spectral 
function. As a first approximation we could pick only the diffusion mode. However, the 
spectral function and the residue have the same sign, so the approximation would clearly 
fail at the first zero. This means that in order to have a sensible approximation for larger 
values of the momentum, we must take into account higher modes. We will see that this 
solves the problem. 

Following the discussion in section ||, we can compare the results from the exact re- 
tarded density-density correlator Gtt with the approximation taking the first four quasinor- 
mal modes or just the hydrodynamic mode. In order to fix the analytic term C4 in (|3.7|). 
we equate the values of the exact and the approximate Green functions at zero frequency 
and q = 0.2, giving C4 ~ 3.9. In figure [l] we can see that the approximation for fixed 
momentum is very good in the interval to < 1, even for larger values of the momentum. 




Figure 1: Imaginary (left) and real (right) parts of the retarded Gtt correlator as a function of 
the frequency at q = 0.2 (black) and q = 0.4 (grey). The dotted line is the hydrodynamic mode 
contribution, the solid line is the exact solution and the dashed line is the four-mode approximation. 

We also present the residues for the complex momentum modes, in figures [| and |3|. 
Again these are the residues of the scalar functions ILt L- We find that contrary to the 
quasinormal modes, the hydrodynamic mode does not decay at high frequencies (see figure 
|3|). In [12] the locations of the complex momentum modes in the R-charge diffusion channel 
have been studied. It was found numerically that the real part of the complex momentum 
approaches ui whereas the imaginary part becomes smaller and smaller at high frequency 
for all the modes. Furthermore we argued in the introduction that the high frequency and 
momentum limit can be exchanged with the zero temperature limit T —* 0. The theory 
at T = recovers Lorentz symmetry and in the case at hand even conformal symmetry. 
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This restricts all signal propagation at T = automatically to the light cone, showing that 
indeed vp = 1. 



3 : 
« 2 - 




1 2 3 4 5 6 




2 ^F^-n^i-r-T--rT'-rr- ; "T"-; 

1 2 3 4 5 6 



Figure 2: (Left) Real and imaginary parts of the residues for the first four complex momen- 
tum modes in the transverse component IIt- (Right) Idem for longitudinal component IIl. The 
numerical values are normalized by N 2 T 2 /8 and the square of the mode number. 
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Figure 3: Real and imaginary parts of the residues for the diffusion complex momentum mode. 
The numerical values are again expressed in terms of N 2 T 2 /8. 



5. Energy-momentum tensor correlators 

The energy-momentum tensor is a conserved quantity 

= , (5.1) 

we can identify E = J d 3 xToo and Pi = J d 3 xToj as the conserved energy and momentum. 
Since they are conserved and cannot be dissipated away, they will slowly spread through 
the plasma or they will be displaced between different regions. This is described by the hy- 
drodynamic shear and sound modes. The response to an external perturbation represented 
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by the source j^ u (t,x) is given by 

(6T^(t,x)) = - J diMVG^^i-^x-xOrV,*') , (5-2) 

where Gr is the retarded two-point function. For low frequencies and momentum (large 
times and distances), the response is dominated by the hydrodynamic modes. 

The energy-momentum tensor T^ u in the gauge theory is dual to the five-dimensional 
metric qmn- We work in the linear approximation, so we will consider small fluctuations 
h-MN on top of the background metric of a planar black hole in AdSs space ( |1 . 1| ) . We fix 
the gauge to h r M = and expand the metric components in plane waves hu U (t,yL,r) = 
e —tut-Hq?cf l f r y j£ j s convenient to use quantities that are gauge-invariant under the 
residual diffeomorphism transformations and group them according to their representation 
under the rotational group. In this way, we can distinguish three different components, 
scalar, vector and tensor, that correspond to the sound, shear and scalar channels in the 
dual theory [41]. Assuming that the momentum is along the z direction, and using the 
dimensionless frequency and momentum 27rT(tD, q) := (uj, q), the gauge- invariant quantities 
are 



2 



(Shear) Z u = c\h ti /r 2 + tvh zi /r 2 , i = x,y , 

(Sound) Z 2 = q 2 h tt /r 2 + 2tt>qh tz /r 2 + Vo 2 h zz /r 2 + (q 2 - tt> 2 + rf'(r)/2)(h xx + h yy )/r 
(Scalar) Z 3 = h xy /r 2 . (5.3) 

The retarded Green function of the dual gauge theory is computed as the boundary action 
of classical solutions satisfying infalling boundary conditions at the horizon, taking into 
account the subtleties of the Lorentzian formulation [7,42]. The poles of the retarded 
correlator correspond to the quasinormal modes, that are normalizable solutions at the 
AdS boundary. The method to find the quasinormal modes is explained in appendix [A|. 
The Green functions G\, G 2 , G3 found using the Z\,Z 2 , Z3 components are scalar quantities. 
They are the coefficients of the tensor projectors depending on uj and q in which energy- 
momentum correlators can be decomposed once Lorentz invariance is reduced to rotational 
invariance (c.f. [41]). For instance, 



Gtx.tx — tt - T~ — 9 G\ , (5-4) 



G t t,tt — - , 2 q — JTTTT G 2 , (5.5) 
3 (or — q A ) z 

G X y,xy = — G3 . (5-6) 

We also compute the residues of the different modes. Notice that we are using the functions 
G\, G 2 and G3, so in order to recover the right residue for the different components of the 
energy-momentum tensor correlators, one should take into account the appropriate factors 
coming from the tensor projectors. 

In the vector and scalar channels, the lowest quasinormal mode shows the right be- 
haviour at low momentum to be identified with the shear (uj = —iDq 2 ) and sound (uj = 
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v s q — iT s q 2 ) modes respectively. In figure f| we compare our numerical results with the 
second order hydrodynamic approximations of [23]. From the Kubo formula of the tensor 
component G xyjXy the second order corrections including a term r] rn oj 2 have been identified 
with 

2~log2 

m = • (5.7) 

According to the hydrodynamic expansion, the first corrections to the shear and sound 
poles would be 

■ V 2 • f V \ 2 4 
Wshear = ~^<? ~ 1 [^f ) T n 1 , 

^sound = v s q - iT s q 2 + ^ \ v 2 s r n - y^j <? 3 . (5.8) 

With rj the shear viscosity, s the entropy density and v s and T s the sound velocity and the 
sound attenuation. However, in the holographic computation the dispersion relation for 
the shear mode differs from the second order hydrodynamic result. This discrepancy can 
be attributed to the necessity of going to third order hydrodynamics to capture the right 
0(q 4 ) correction to the shear pole. In terms of the dimensionless variables the holographic 
theory gives the next to leading order corrections in the dispersion relations 

.q 2 , (l-log2) 4 

tfshear = ~i — ~ « T q , 

q .q 2 3-21og2 , . . 

K — = vi-'i + -wr' ■ (5 - 9) 

When we compare the numerical results with these (fig. ^), we find a very good agreement 
U P to q < 1. Similar comparisons between the exact values and first order hydrodynamics 
or the first corrections can be found in [11,23,41]. It seems justified then to consider just 
the first hydrodynamic corrections to describe the properties of the shear and sound modes 
up to frequencies of the order of the temperature u, q ~ T. Note however, that the q 
term that arises from second order hydrodynamics obviously does not approximate the 
numerical result very well, clearly indicating that second order hydrodynamics is not valid 
to predict the q terms as already pointed out in [23]. 

We have also computed the complex momentum modes that describe the penetration 
depth of a perturbation inside the plasma as a function of the frequency [12]. Here we 
show the complex momentum modes for the sound channel in figure || while for the shear 
channel can be found in reference above. The results for the residues are displayed in figures 
|and|. The zero frequency value corresponds to the inverse of the screening masses in the 
plasma. Again we find that the imaginary parts become smaller with increasing frequency 
and that the real parts approach the light cone u/q — ► 1 indicating a front velocity vp = 1 
as in the other channels (see Appendix ^]) . 

5.1 Residues and higher thermal resonances 

In figures ^ and || we plot the residues of the lowest four quasinormal modes in the shear and 
sound channels. The residue in the shear mode (figure ^) shows a decaying and oscillatory 
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Figure 4: (Left) The numerical value of the shear mode [solid black] compared with the correct 
expression [dotted gray] and the value from second order hydrodynamics [dashed gray]. (Middle) 
The numerical value of the real part of the sound mode [solid black] compared with the second 
order hydrodynamic approximation [dotted gray]. (Right) The numerical value of the imaginary 
part of the sound mode [solid black] compared with the second order hydrodynamic approximation 
[dotted gray]. 
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Figure 5: (Left) Real and imaginary parts of the first four complex momentum modes in the sound 
channel. (Right) Idem for sound mode. 



pattern with increasing momentum, similar to the diffusion mode of the R-current. We can 
apply the same arguments here, at short enough wavelengths the first quasinormal mode 
will dominate the late time behaviour and the hydrodynamic approximation will not be 
valid. The sound mode on the other hand has a non-vanishing residue at high momentum, 
and its pole moves away from the real axis only for low values of the momentum, but then 
it behaves as an ordinary quasinormal mode, so it always dominates the late time response 
of the system. 

Using the numerical values of the residues, we can study the hydrodynamic time ( |2.6| ) 
and length fl2,10| ) scales for the shear and sound modes. Notice that the residues of the 

(2) 

sound channel R n (q) we have computed are evaluated at the quasinormal frequencies and 
correspond to the 'scalar' Green function G*2. From fl5.4|), we can see that in order to 
expand Gtt, tt as a sum over quasinormal modes we should use the residues 



R (tt,tt) _ ± 



3 (w£(qr) 



.,2^2 



(5.10) 
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Figure 6: (Left) Real and imaginary parts of the residues for the first four complex momentum 
modes in the shear channel. (Right) Idem for sound channel. 




ro m 



Figure 7: (Left) Real and imaginary parts of the residue for the shear mode. (Right) Idem for the 
sound mode. 



In the shear channel we are computing the residues (q) of the function G\ evaluated 
at the quasinormal frequencies. Using (5A), the residues of Gt x ,tx should be 

R^ tx) = \ fX 2 4 1} (g) • (5.H) 

However, contrary to the cases of the density-density correlator of the R-current or the 
Gtt,u correlator of the sound channel, the large frequency behaviour does not asymptote 
to a momentum-dependent constant (~ q 2 or q 4 ), but to ~ q 2 (uo 2 - q 2 ). This means 
that there should be an implicit uj 2 dependence in the residues, but this is not shown by 
our computation. A way to partially recover the right dependence is to use the following 
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definition 

R ( t' tx) = l ( 2 fr g ' 2 ci 2 R^{ q ). (5.12) 

We consider the same kind of perturbations as for the R-current diffusion mode in 
section |3], one is a spatial plane wave of momentum q that is switched on only during a 
finite amount of time At. It sources the transverse and longitudinal momentum densities, 
whose response is given by Gt x ,tx and Gu,tt- The other perturbation is a pulse of size Az 
localized in the z direction that oscillates with frequency uj. We study the response of the 
transverse and longitudinal momentum currents given by for G ZXtZX and G zzzz . The results 
for the shear mode are qualitatively the same as for the diffusion mode. For a short-lived 
time-localized source there is a minimal hydrodynamic time 27tTth ~ 1.34. The maximal 
value of the momentum beyond which the relaxation to equilibrium is dominated by higher 
modes at late times is q ~ 2.6ttT. For a space-localized source, for large enough Az > 1/T 
and small enough frequencies uj < 0.6ttT, the screening is well described by the complex 
momentum shear mode. For smaller sizes or higher frequencies in general the shear mode 
dominates only at a finite distance from the source in > 0, unless the frequency is very 
low. 

The sound mode behaves qualitatively different because its residue does not vanish 
at any value of the momentum. It always dominates the relaxation to equilibrium for 
long enough times or large enough distances as far as our computation can show. Higher 
modes dominate the response for a finite time after switching off the source for momenta 
q > 1.27rT. There is no minimal time for the sound mode, but since for large frequencies 
the sound mode behaves as any other mode, the behaviour will not be hydrodynamic until 
some finite time has passed if the source is very short-lived. If the source is localized in 
space, higher modes dominate at a finite distance for values of the frequency to > 2.5ttT. 

We can also study the validity of the hydrodynamic regime through the spectral func- 
tion. As we found for R-current diffusion, the first zero of the shear residue implies a 
change of sign for the approximation with the shear mode alone, hence a failure of the hy- 
drodynamic approximation. This happens around q ~ 1.3. We find that also the spectral 
function restricted to the sound mode alone changes sign but it is not related to a 'cross- 
ing'. It coincides approximately with the change of behaviour of the residues that can be 
observed in figure ^, close to q ~ 1.1. Therefore, in order to have a consistent description 
of the system for larger values of momentum, it is necessary to take into account higher 
modes. 

Using these results for the residues and the frequencies of the first quasinormal modes 
we can build an approximation to the retarded correlator as explained in section ||. For 
the sound channel, the properly defined Green function actually includes the boundary 
terms found in [41] added to the Frobenius series approximation. Otherwise, an unphysical 
singularity would appear at uj = ±q. We examine the Gtt,tt Green function and find 
that the approximation and the exact function differ by a real constant C ~ 0.374. Once 
this term is included, they agree for a fair interval of frequencies and momenta (fig. 10). 
Even for values of u, q smaller but comparable to the temperature, the hydrodynamic 
mode alone gives a remarkable good approximation for the spectral function. We can 
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Figure 8: (Left) Real and imaginary parts of the residues for the first four quasinormal modes in 
the shear channel. (Right) Idem for the sound channel. The numerical values are normalized by 
7r 2 7V 2 T 4 and the square of the mode number. 




Figure 9: (Left) When the shear mode hits a special value, there is a qualitative change in the 
behaviour of the quasinormal modes. Zeroes of the residue of the shear mode roughly coincide with 
the crossing of the shear mode and the quasinormal modes. (Right) The residue of the sound mode 
behaves as the rest of the quasinormal modes. The numerical values are normalized by n 2 N 2 T A 



see that for q = 0.2, the hydrodynamic approximation and the exact result are virtually 
indistinguishable for to < 1. 

Other correlators related to the sound channel can be found from Gtt,tt using the 
explicit expressions for the tensor projectors [41]. 



We also examine the Gt x ,tx component of the shear channel, the results are in fig. 11 
We have to add an analytic piece ~ 0.9q 2 ro to the spectral function. To the real part of 
the correlator we need to add a more complicated term ~ 1.89q 4 - 0.16q 2 - 1.77q 2 tt> 2 . The 
plots made in figure 11 have been done using the definition ( |5.12j ), we observe very good 
agreement for to < 1, q < 0.6. 
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Figure 10: Imaginary (left) and real (right) parts of the retarded G u , tt correlator as a function 
of the frequency at q = 0.2 (black) and q = 0.4 (grey) for the real part and q = 1 (grey) for the 
imaginary part. The dotted line is the hydrodynamic mode contribution, the solid line is the exact 
solution and the dashed line is the four-mode approximation. 
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Figure 11: Imaginary (left) and real (right) parts of the retarded Gt x ,tx correlator as a function 
of the frequency at q = 0.2 (black) and q = 0.4 (grey). The dotted line is the hydrodynamic mode 
contribution, the solid line is the exact solution and the dashed line is the four-mode approximation. 

6. Discussion 

We have studied hydrodynamics in the strongly coupled N = 4 gauge theory based on 
linear response theory using the AdS/CFT correspondence. As emphasized we understand 
hydrodynamics here as the effective theory resulting from integrating out the higher quasi- 
normal modes and keeping only the contribution of the hydrodynamic modes, i.e. those 
modes whose quasinormal frequency vanishes at zero momentum. We also recall that we 
can consider this to be a definition of the system being in local thermal equilibrium. One 
of the important findings is that the hydrodynamic approximation defined in that way has 
its breakdown built into it: we saw that the diffusion mode in the shear channel crosses 
the lowest quasinormal mode at around q ~ 1.3. From that value on it is the lowest non 
hydrodynamic mode that determines the late time behaviour. In fact, already slightly 
before the spectral function of the hydrodynamic mode ceases to be positive. We interpret 
this as a signal for the breakdown of the effective theory based on the hydrodynamic mode 
alone. Similarly, in the sound channel we found that the spectral function of the sound 
mode switches sign at q ~ 1.1 and again we take this as a breakdown of a putative effective 
theory based on the sound channel alone. The full spectral functions (or even the approx- 
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imation keeping the contributions from only a few quasinormal modes) behave perfectly 
reasonable at these points. 

An important role in our investigations have been played by the residues. We saw 
that the shear diffusion mode and the R-charge diffusion are very similar. Both decay 
in an oscillatory pattern with decreasing frequency and numerically we find that they 
decouple for momenta q > 1. As we explain in appendix B this can be understood from 
the analyzing the in-falling boundary conditions at the special values of the frequency 
ro = —in with n £ N. That both diffusion modes show such similarities leads one naturally 
to suspect that this might be a universal property of holographic hydrodynamic diffusion 
modes. At the moment it is not clear how such a conjectured universal behaviour could 
be proved. However it seems doable and an interesting task to generalize the calculations 
of the residues presented here and in [29] to other holographic gauge theories and check if 
the diffusive modes behave in the same way [43]. A clue of why this might happen comes 
from causality: in order to preserve causality and at the same time reproduce (first order) 
hydrodynamics in the long wavelength limit some drastic modification at short wavelengths 
is definitely necessary. 

As we have seen, the complex momentum modes needed for the calculation of the front 
velocity do behave perfectly causal and their residues do not decouple at high frequencies. 
Although in the small frequency/momentum limit the lowest complex momentum mode 
can be computed from the analytic continuation of the quasinormal hydrodynamic mode 
ro = — iq 2 this is not so for short wavelengths and high frequencies. This does not come as 
a surprise, since only rotational invariance is preserved in the finite temperature theory, so 
non-analyticities between oj and q dependence are expected, like uj/q. An example of this is 
the fact that the two limits u> — > and q — > of the retarded correlator do not commute in 
general. Also, in order to do the analytic continuation of the Green function properly, all 
modes must be taken into account, only in the u — > 0, q — > limit the hydrodynamic mode 
dominates. It is interesting to consider the hydrodynamic time scale we found in the energy 
density correlator, 2ttTtu ~ 1.34. At RHIC temperatures T ~ 300 MeV this translates into 
a very short time th = 0.14fm/c. 3 One might take this as an indication for an extremely 
fast thermalization time. Of course, thermalization at RHIC includes processes that are 
outside the regime of linear response considered here, so the short hydrodynamic scale can 
at best describe a late stage of thermalization. 

Another important point was to see how well the retarded Green's functions can be 
approximated by keeping only a few quasinormal modes. We found that analytic pieces 
related to the non-convergence of the sum over quasinormal modes played an important 
role. Recently, attempts of reconstructing the quasinormal mode spectrum from (much 
easier to compute) holographic spectral functions have been made in [44]. We think that 
our observations might also be useful to gain a better control on such procedures. 



3 For the R-charge density this is even around four times shorter. This time has been erroneously reported 
to be 0.3 fm/c in [29]. 
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A. Method 

A.l Energy- momentum tensor 

In the following we will use the coordinate u = r^/r 2 such that the horizon sits at u = 1 
and the boundary at u = 0. The equations of motion for the diffeomorphism-invariant 



quantities defined in (^3) are [41] 



(tr 2 - q 2 f(u))f{u) - utv 2 f'(u) tu 2 - tffju) 

Z\ H 77 w 277 a — ZjT\ Z i H TTvi — Zl = > (A.la) 

uf\u)(q 2 f(u) - tt) 2 ) uf(u) 2 

_ 3ro 2 (l + ^ 2 ) + q 2 (2n 2 -3n 4 -3) 

2 u/(u)(3tt> 2 + q 2 (u 2 -3) 2 

3ro 4 + q 4 (3-4u 2 + n 4 ) + q 2 (4n 5 -4n 3 + 4ttJ 2 n 2 -6lD 2 ) 

^ 7? \2^o^2 — 27 _ 2 ovi ^2 = 0, (A.lb) 

uj(uy{3tt> z + q^(u z - 3)) 

„ 1 + u 2 , to - q 2 /(^) 
Z 3+— ??-T Z 3 + > / v 2 Z 3 = • (A.lc 

Both the boundary and the horizon are regular singular points. At the horizon (u = 1), 
there are two possible local solutions 

Z (a) ~ (1 - n)~-/ 2 ^ + (1 - u)^ 2 <fffi , a = 1, 2, 3. (A.2) 

In order to compute the retarded Green function, we must pick infalling boundary condi- 
tions (p°^ = 0. According to the holographic dictionary the retarded Green function can 
be computed as the ratio of the connection coefficients that relate the local solution at the 
horizon with ingoing boundary conditions to the non-normalizable (A a ) and normalizable 
(B a ) solutions at the boundary (u = 0) 

The Green functions in the different channels are determined by three scalar functions G( a ) 
given by the ratio of the connection coefficients 

G (a) = -n 2 N 2 T^. (A.4) 
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The quasinormal modes are normalizable solutions where At a ) = 0. The ratio Bu\/A( a ) 
follows from demanding that the solution is smooth at a matching point in the interior of 
the interval (0, 1) 



B 



(«) 



A, 



Z (a) ( Z (a) )' ~ ( Z (a) ) "(ft) 

We have computed the Frobenius series up to order 50. Matching the series expansions, 



(A.5) 



we see that the ratio (A.5) remains constant for a fair interval in the radial coordinate. We 
have chosen x = 0.53 to evaluate the ratio. 

The residue R n for the quasinormal mode to n can be computed as 



Rn 



JL 






dm 


Of) 


tv=m n _ 



(A.6) 



A. 2 Global current 



The equations of motion for the gauge invariant combinations (section ^ of the plane- wave 
vector field perturbations are [41] 



f\u) , to' 



f(uW 



(uf(u) 



/(«)(t» 2 -/(«)q 2 



E'r + 



to' 



f(uW 



uf(u) 



E T = , 



E L = 



(A.7a) 
(A.7b) 



Defining (a) := (T, L) as the two gauge-invariant components components, we can fol- 
low the same procedure as with the energy-momentum tensor components. The infalling 
solution can be expanded at the boundary in the non-normalizable and normalizable modes 



£&)(«) = Aa) E} a) (u) + B (a) Ef a) (u) 



(a) 



(Ai 



where E^(u) ~ u. The retarded Green function is determined by the longitudinal and 
transverse polarization 



8 A( a ) 

which are proportional to the ratio between the connection coefficients 



(A.9) 



B 



(a) 



Z?AE}J'-{E?JE} 



A 



(a) 



(A.10) 



We have computed the Frobenius series up to order 50. Matching the series expansions, 
we see that the ratio ( [A. 10 ) remains constant for a fair interval in the radial coordinate. 
We have chosen x = 0.53 to evaluate the ratio and have checked that the spectral function 
agrees with previous numerical (for non-zero momentum) and exact (for zero momentum) 
results [36,45,46]. The quasinormal modes correspond to the frequencies where there is a 
pole .4 = 0. We can apply equation (|A.6|) to compute the residues. 
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B. Zeroes of hydrodynamic residues 

We have seen that the residues of the diffusion and the shear mode have an oscillatory 
behaviour with the momentum. In this section we show how to find the location of the 
zeroes of the residues. We will use that the equations of motion for the for the vector field 
and vector component of the metric are Heun equations (c.f. [12, 33]) 4 

»"« + (i + + ^2) + ^W* -2) »w = • < R1 » 

In the case of the vector field, we define Vq(z) = A'q(z), Vl(z) = A' L (z). In the 
x = 1 — z 2 coordinate, the Heun equations are found using the new variables 



V (x) = x~ ixo ' 2 {x - l)^ 2 (x - 2)- m / 2 y{x) , 
V L (x) = x^-^ix - l)^ 2 (x - 2)- 1 - K, l 2 y{x) . 
In both cases we find the same parameters for the Heun equation 

a = -£(l + t), / 9 = 2-(|(l + i)) , Q = q 2 -(l + 3z)|-(2-z)^ 
7 = 1 — iro, 8=1, e=l-tn. 



(B.2) 



(B.3) 



For the shear component we use the gauge-invariant variable ipv proposed in [47]. The 
Heun equation is found for the new variable 

Vv(x) = x- iro/2 (x - l) 3/4 (x - 2)- ro / 2 y(x) , (B.4) 

with parameters 

aP=^(l + i)(^(l + i)-3) , Q = q 2 -(l + 5i)^-(2-i)^- , 

7 = 1 — itt> , 8 = 2, e = 1 - to . 
The coefficients of the Frobenius series at x = should satisfy the recursion relation 

2(n + 2)(n+ 1 + 7K+2 + A n (to , <\) a n+1 +B n (tv,q)a n = , n>0, (B.6) 

where 

A n (to, q) = -((n + 1)(2J + e + 3(n + 7)) + Q) , (B.7) 
B n (to,q) = (n + a){n + P) , (B.8) 

and 2701 — Qao = 0. This recursion relation has in general a unique solution. However, 
when the continued fraction 

a n A n (u,q) + r n+ i 



4 There is a factor of two difference with the conventions used here for the frequency and the momentum. 
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converges, Pincherle's theorem states that an extra solution to the three term recursion 
relation ( |B,6D exists [9,10,48]. This is equivalent to finding a solution that is analytic 
both at the boundary and at the horizon. Notice that at the horizon (x = 0), the critical 
exponents are and 1 — 7 = iro, so the analytic solution usually corresponds to the 
first solution, that is the one associated to infalling boundary conditions. However, when 
ro = — ik, with k = 1, 2, . . . , the analytic solution will be in general the outgoing solution, 
since it corresponds to the largest integer critical exponent, while the infalling solution 
will have a logarithmic contribution. This is reflected in the recursion relations, since the 
coefficient of the a n +2 term vanishes when n + 2 = k, implying that the series start at x k . 
The reason for this is that the first k recursion relations form a closed subsystem of linear 
equations for k unknowns. Then, the general solution to the Heun equation will be of the 
form 

00 / 00 00 \ 

y{x) = a 1 x k ^2a n x n + a 2 l^2b n x n + c{q)log{x)x k ^2a n x n ) . (B.10) 
n=0 \n=0 n=0 / 

In principle, there could be special values of the momentum where the coefficient of the 
logarithm vanishes c(q) = and the two solutions at the horizon are analytic. This actually 
happens when the rank of the subsystem of the first k recursion relations is zero. In that 
case, we can find trivially two solutions satisfying the recursion relations without having 
to check the convergence of ( |B.9[ ). 

The diffusion and shear modes are located at negative imaginary values of the fre- 
quency, larger as the momentum is increased. This implies that for some large enough 
value of the momentum, the frequency will have the special value ro = —ik. However, the 
mode has to be analytic at the boundary and the horizon, so the value of the momentum 
when the special value is reached must be determined by the condition c(q) = 0. When we 
compare the analytic result with the numerical computation, we find that this is indeed 
the case, the first points are (ito, q 2 ) = (1, 1/2), (2, - 1), (3, y/E - 3/2) for the diffusion 
mode and (itt>,q 2 ) = (1,VE), (2,3.2266), (3,3.91764) for the shear mode. On the other 
hand, at the special point the analytic solution that is found as the limit of the hydro 
mode ro — * —ik is a linear combination of the normalizable and non-normalizable modes. 
This implies that the pole disappears from the retarded Green function or in other words, 
that the residue is zero. 

C. Front velocity 

Wave propagation in dispersive media has been studied long ago in the classic work of 
Brillouin and Sommerfeld [49-51]. It has been pointed out there that the group velocity 
v g = duj/dq is not a reliable indicator if one wants to study the question of how fast can 
a signal be transmitted through the dispersive medium. In fact it is the so called front 
velocity which limits the speed of propagation of a signal through the medium. The front 
velocity is defined as the velocity with which the onset of a signal travels. In dispersive 
media this is not yet sufficient to guarantee that the signal travels with this speed and 
therefore one also has to define a signal velocity, which is the speed with which practically 
usable signals travel. This signal velocity is always smaller than the front velocity. For 
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matters of principle, i.e. answering the question if causality is preserved it is therefore 
the front velocity that is important. We will briefly review here the reasoning leading to 
the definition of the front velocity. We want to study how fast a perturbation can travel 
through the medium, for this purpose we will switch on a periodic signal with frequency 
v at time t = 0. We model this by a source of the form @(t)e~ tut 5(x). 5 To compute 
the response of the system we expand the retarded Green function in its poles in the 
complexified momentum plane (for simplicity we will also restrict our considerations to 
effectively 1 + 1 dimensions) 

<*(*,*)> = - f to^yRnM) i e - Mqx {C1) 

J 2ir 2tt ^-^ q — q n u — v + ie 

The poles in the momentum k have real and imaginary parts q n {w) = <2^(w) + iq^(u) and 
from [12,33] we know that they have to lie symmetrically in the first and third quadrants. 
We consider therefore only the region with x > since the response in x < is just the 
mirror image. Assuming x > we have 

(*(M)> = - [ ^Vi?„(^n) l —e~ w i+vz^e-*"*. (C.2) 



To pick up a non vanishing signal we have to close the contour of integration now in the 
lower half plane. This demands however that ]im CJ -^ 00 (t — x qn ^ ) > showing that we 
can pick up the first front of the signal only at space time points where x/t < vp, where 
we have defined the front velocity 6 

vp := lim — p . (C3) 

uj^oo q u 

It follows that causality is preserved if the front velocity is smaller that the speed of light 
vf < c. Results for the complex momentum modes for the shear channel and the diffusion 
channel in [12,33] and for the sound channel in figure [|, show that the imaginary parts 
go to zero if we extrapolate to the large frequency limit, i.e. u — * oo, while the real 
part approaches from below the dispersion relation co = q. Moreover, the evolution of the 
hydrodynamic modes in this limit is the same as the evolution of the higher modes, finding 
numerically that lim^^oo — = 1. This provides a numerical proof that the strongly coupled 
M = 4 plasma behaves causally even when we take into account the hydrodynamic modes. 

References 

[1] Juan Martin Maldacena. The large N limit of superconformal field theories and supergravity. 
Adv. Theor. Math. Phys., 2:231-252, 1998. hep-th/9711200. 

[2] S. S. Gubser, Igor R. Klebanov, and Alexander M. Polyakov. Gauge theory correlators from 
non-critical string theory. Phys. Lett, B428:105-114, 1998. hep-th/9802109. 



5 For simplicity we chose to localize the source in space, it is not important for the general argument. 

6 Since this implies a linear relation between to and q for large frequencies and (real) wave numbers it is 
usually written as vf = lim^oo — . It also follows that the diffusion equation violates causality since from 
the dispersion relation we find ui — 2Dq 2 and therefore uj/q = 2Dq with no finite limit for the front velocity. 



- 25 - 



Edward Wittcn. Anti-dc Sitter space and holography. Adv. Theor. Math. Phys., 2:253-291, 
1998. hep-th/9802150. 

S. S. Gubser, Igor R. Klebanov, and A. W. Peet. Entropy and Temperature of Black 
3-Brancs. Phys. Rev., D54:3915-3919, 1996. hep-th/9602135. 

Igor R. Klebanov and A. A. Tseytlin. Entropy of Near-Extremal Black p-brancs. Nucl. 
Phys., B475:164-178, 1996. hep-th/9604089. 

Edward Witten. Anti-de Sitter space, thermal phase transition, and confinement in gauge 
theories. Adv. Theor. Math. Phys., 2:505-532, 1998. hep-th/9803131. 

Dam T. Son and Andrei O. Starinets. Minkowski-space correlators in AdS/CFT 
correspondence: Recipe and applications. JEEP, 09:042, 2002. hep-th/0205051. 

Gary T. Horowitz and Veronika E. Hubeny. Quasinormal modes of AdS black holes and the 
approach to thermal equilibrium. Phys. Rev., D62:024027, 2000. hep-th/9909056. 

Danny Birmingham, Ivo Sachs, and Sergey N. Solodukhin. Conformal held theory 
interpretation of black hole quasi-normal modes. Phys. Rev. Lett., 88:151301, 2002. 
hep-th/0112055. 

Andrei O. Starinets. Quasinormal modes of near extremal black branes. Phys. Rev., 
D66:124013, 2002. hep-th/0207133. 

Alvaro Nunez and Andrei O. Starinets. AdS/CFT correspondence, quasinormal modes, and 
thermal correlators in N = 4 SYM. Phys. Rev., D67:124013, 2003. hep-th/0302026. 

Irene Amado, Carlos Hoyos, Karl Landsteiner, and Sergio Montero. Absorption Lengths in 
the Holographic Plasma. JEEP, 09:057, 2007. arXiv : 0706 . 2750 [hep-th] . 

Giuseppe Policastro, Dam T. Son, and Andrei O. Starinets. From AdS/CFT correspondence 
to hydrodynamics. JEEP, 09:043, 2002. hep-th/0205052. 

Giuseppe Policastro, Dam T. Son, and Andrei O. Starinets. From AdS/CFT correspondence 
to hydrodynamics. II: Sound waves. JEEP, 12:054, 2002. hep-th/0210220. 

G. Policastro, D. T. Son, and A. O. Starinets. The shear viscosity of strongly coupled N = 4 
supersymmctric Yang-Mills plasma. Phys. Rev. Lett., 87:081601, 2001. hep-th/0104066. 

P. Kovtun, D. T. Son, and A. O. Starinets. Viscosity in strongly interacting quantum field 
theories from black hole physics. Phys. Rev. Lett., 94:111601, 2005. hep-th/0405231. 

Bcrndt Muller and James L. Nagle. Results from the Relativistic Heavy Ion Collider. Ann. 
Rev. Nucl. Part. Set., 56:93-135, 2006. nucl-th/0602029. 

W. Israel. Nonstationary irreversible thermodynamics: A Causal relativistic theory. Ann. 
Phys., 100:310-331, 1976. 

I. Muller. Zum Paradoxon der Warmeleitungstheorie. Zeitschrift fur Physik A Eadrons and 
Nuclei, 198:329, 1967. 

Stewart J.M. On Transient Relativistic Thermodynamics and Kinetic Theory. Proc. Roy. 
Soc, 357:59, 1977. 

W. Israel and J. M. Stewart. Transient relativistic thermodynamics and kinetic theory. Ann. 
Phys., 118:341-372, 1979. 



- 26 - 



[22] Michal P. Heller and Romuald A. Janik. Viscous hydrodynamics relaxation time from 
AdS/CFT. Phys. Rev., D76:025027, 2007. 

[23] Rudolf Baier, Paul Romatschke, Dam Thanh Son, Andrei O. Starinets, and Mikhail A. 

Stephanov. Relativistic viscous hydrodynamics, conformal invariance, and holography. 2007. 
0712.2451. 

[24] Sayantani Bhattacharyya, Veronika E Hubeny, Shiraz Minwalla, and Mukund Rangamani. 
Nonlinear Fluid Dynamics from Gravity. JHEP, 02:045, 2008. 

[25] Makoto Natsuume and Takashi Okamura. Causal hydrodynamics of gauge theory plasmas 
from AdS/CFT duality. Phys. Rev., D77:066014, 2008. 0712.2916. 

[26] Michel Le Bellac. Thermal Field Theory. Cambridge University Press, 2000. 

[27] Ulrike Kraemmer and Anton Rebhan. Advances in perturbative thermal held theory. Rept. 
Prog. Phys., 67:351, 2004. hep-ph/0310337. 

[28] Jean-Paul Blaizot and Edmond Iancu. The quark-gluon plasma: Collective dynamics and 
hard thermal loops. Phys. Rept, 359:355-528, 2002. hep-ph/0101103. 

[29] Irene Amado, Carlos Hoyos-Badajoz, Karl Landsteiner, and Sergio Montero. Residues of 

Correlators in the Strongly Coupled N=4 Plasma. Phys. Rev., D77:065004, 2008. 0710.4458. 

[30] Joshua J. Friess, Steven S. Gubser, Georgios Michalogiorgakis, and Silviu S. Pufu. Expanding 
plasmas and quasinormal modes of anti-de Sitter black holes. JHEP, 04:080, 2007. 

[31] Suphot Musiri, Scott Ness, and George Siopsis. Perturbative calculation of quasi-normal 
modes of AdS Schwarzschild black holes. Phys. Rev., D73:064001, 2006. 

[32] George Siopsis. Large mass expansion of quasi-normal modes in AdS(5). Phys. Lett., 
B590:105-113, 2004. 

[33] Carlos Hoyos, Karl Landsteiner, and Sergio Montero. Holographic Meson Melting. JHEP, 
04:031, 2007. hep-th/0612169. 

[34] R. A. Konoplya. On quasinormal modes of small Schwarzschild- anti-de-Sitter black hole. 
Phys. Rev., D66:044009, 2002. hep-th/0205142. 

[35] Carlos Hoyos-Badajoz, Karl Landsteiner, and Sergio Montero. Quasinormal modes and 
meson decay rates. Fortsch. Phys., 55:760-764, 2007. 

[36] Robert C. Myers, Andrei O. Starinets, and Rowan M. Thomson. Holographic spectral 
functions and diffusion constants for fundamental matter. JHEP, 11:091, 2007. 
arXiv: 0706. 0162 [hep-th] . 

[37] Angel Parcdes, Kasper Peeters, and Marija Zamaklar. Mesons versus quasi-normal modes: 
undercooling and overheating. JHEP, 05:027, 2008. 

[38] Edward V. Shuryak and Ismail Zahed. Rethinking the properties of the quark gluon plasma 
at T approx. T(c). Phys. Rev., C70:021901, 2004. hep-ph/0307267. 

[39] Gerald E. Brown, Chang-Hwan Lee, Mannque Rho, and Edward Shuryak. The anti-q q 
bound states and instanton molecules at t >= T(C). Nucl. Phys., A740:171-194, 2004. 
hep-ph/0312175. 

[40] Edward V. Shuryak and Ismail Zahed. Towards a theory of binary bound states in the quark 
gluon plasma. Phys. Rev., D7O:054507, 2004. hep-ph/0403127. 



-27- 



[41] Pavel K. Kovtun and Andrei O. Starinets. Quasinormal modes and holography. Phys. Rev., 
D72:086009, 2005. hep-th/0506184. 

[42] C. P. Herzog and D. T. Son. Schwinger-Keldysh propagators from AdS/CFT correspondence. 
JHEP, 03:046, 2003. hep-th/0212072. 

[43] Work in progress. 

[44] Robert C. Myers and Aninda Sinha. The fast life of holographic mesons. 2008. 0804.2168. 

[45] Derek Teaney. Finite temperature spectral densities of momentum and R-charge correlators 
inN = 4 Yang Mills theory. Phys. Rev., D74:045025, 2006. hep-ph/0602044. 

[46] Pavel Kovtun and Andrei Starinets. Thermal spectral functions of strongly coupled N = 4 
supersymmetric Yang-Mills theory. Phys. Rev. Lett., 96:131601, 2006. hep-th/0602059. 

[47] Hideo Kodama and Akihiro Ishibashi. A master equation for gravitational perturbations of 
maximally symmetric black holes in higher dimensions. Prog. Theor. Phys., 110:701-722, 
2003. hep-th/0305147. 

[48] E. W. Leaver. An Analytic representation for the quasi normal modes of Kerr black holes. 
Proc. Roy. Soc. Lond., A402:285-298, 1985. 

[49] A. Sommerfeld. Ein Einwand gegen die Relativtheorie der Elektrodynamik und seine 
Beseitigung. Physikalische Zeitschrift, 8:841, 1907. 

[50] L. Brillouin. Uber die FortpBanzung des Lichts in dispergierenden Medien. Annalen Phys., 
10:203, 1914. 

[51] L. Brillouin. Wave Propagation and Group Velocity. Academic Press, 1960. 



- 28 - 



